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Within the scope of the locally rotationally symmetric (LRS) Bianchi type-I cosmological 
model the role of spinor field on the evolution of the Universe is investigated. In doing so, we 
have considered a polynomial type of nonlinearity. It is found that, depending on the sign 
of the self-coupling constant, the model allows either an accelerated mode of expansion 
or an oscillatory mode of evolution. While the non-diagonal components of the energy- 
momentum tensor of the spinor field in the case of a full Bianchi type-I model lead to the 
vanishing mass and nonlinear term in the spinor field Lagrangian, in the case of an LRS 
Bianchi type-I model neither the mass term nor the nonlinear term of the spinor field vanish. 
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I. INTRODUCTION 


Recently, after some remarkable works by different authors IH15IL showing the important role 


that spinor fields play on the evolution of the Universe, the situation began to change. This change 
of attitude is directly related to some fundamental questions of modern cosmology: (i) problem of 
initial singularity; (ii) problem of isotropization and (iii) late time acceleration of the Universe. 

Given the role that the spinor field can play in the evolution of the Universe, the question that 
naturally emerges is, if the spinor field can redraw the picture of evolution caused by perfect fluid 
and dark energy, is it possible to simulate perfect fluid and dark energy by means of a spinor field? 
An affirmative answer to this question was given in a number of papers [16h2u]. In those papers, a 
spinor description of matter, such as a perfect fluid and dark energy, was given and the evolution of 
the Universe, given by different Bianchi models, was thoroughly studied. In almost all the papers 
the spinor field was considered to be a time-dependent function and its energy-momentum tensor 
was given by the diagonal elements only. 

Some latest studies show that because of the specific connection with the gravitational field 
the energy-momentum tensor of the spinor field possesses non-trivial non-diagonal components as 
well, and these non-zero non-diagonal components of the energy-momentum tensor play decisive 
a role in the character of the geometry of space-time as well as on the components of the spinor 


field |21|-!24|. 


In this paper we study the evolution of the Universe filled with spinor field within the scope 
of a locally rotationally symmetric (LRS) Bianchi type-I (BI) cosmological model. It should be 
noted that ageneral BI model in the presence of a nonlinear spinor field duly evolves into a LRS 


BI model [22]. In this paper we thoroughly study the role of a spinor field in the evolution of the 


Universe given by a LRS BI model. Here we also consider a more general type of nonlinearity. 


II. BASIC EQUATIONS 

In this paper we plan to study the evolution of the Universe given by a LRS BI anisotropic 
cosmological model filled with nonlinear spinor field. 

The LRS BI model is the ordinary BI model with two of the three metric functions being equal 
to each other and can be given by 


ds — dt —ci i \dx dx J — ci^dz ■ 


(2.1) 
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with a i and a 3 being functions of time only. 

The nontrivial components of the Einstein tensor corresponding to metric (12.11) are 
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( 2 . 2 a) 

( 2 . 2 b) 

(2.2c) 


Keeping in mind the symmetry between y and y we choose the symmetrized Lagrangian [250 
for the spinor field as L5]: 


L = 


VfVuV-VnWfV 


-m sp yy-F, 


(2.3) 


where the nonlinear term F describes the self-interaction of a spinor field and can be presented as 
some arbitrary function of invariants generated from the real bilinear forms of a spinor field. We 
consider F = F(K), with K taking one of the following expressions {/, J, I + J, I — J}. It can be 
shown that such a choice describes the nonlinearity in its most general form. 

Here Vp is the covariant derivative of spinor field 


dy 


dy 




(2.4) 


with F u being the spinor affine connection. In (12.31) y’s are the Dirac matrices in curve space-time 
and obey the following algebra: 

/y + r v y M = 2^ v ( 2.5) 


and are connected with the flat space-time Dirac matrices y in the following way 

8fiv(x) = e a Jx)e h v (x)ri ah , y^x) = e a Jx)%, 


( 2 . 6 ) 


where ry,/, = diag — (1, — 1, — 1, — 1) and e a ^ is a set of tetrad 4-vectors. The spinor affine connection 
matrices rip (x) are uniquely determined up to an additive multiple of the unit matrix by 


^ nJv — -=r-x ~^%Yp — r^yv + y v rp — 0, 


(2.7) 


with the solution 


- v 1 a 


V n 
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The spin affine connection corresponding to LRS BI metric (12.11) can be written explicitly as 

r 0 = o, r 1 = ^ff, r 2 = r 3 = ^ff, (2.9) 

Varying (12.31) with respect to iff (iff) one finds the spinor field equations: 

¥ - m sp ¥ - 2 F k (SKj + iPKjf)¥ = 0, (2.10a) 

+m sp \fr + 2F K \ff(SK 1 + iPK J y 5 ) =0. (2.10b) 

Here we denote Fk = dF/dK, Ki = dK/dl, and Kj — dK/dJ. 

The energy-momentum tensor of the spinor field is given by 

Tjl = ^g pyr UfYn^v¥+¥Y^n¥-^n¥Yv¥-^v¥Yn¥^ ~^L sp ( 2 . 11 ) 

where L sp in view of ( 12 .lOab and ( 12.1 Obi) can be rewritten as 

L, v = l -[¥Y^n¥-^nWY p ¥] -m sv ¥¥-F(K) 

= ^ VifVn ¥ - >n sp ¥} ~ \ [V M ¥Y P + ™s P ijf\v-F{K), 

= 2 F k ( IK / + JKj) —F — 2 KF k -F(K). (2.12) 

Further inserting (12.41) into (12.1 II) the energy-momentum tensor of the spinor field can be written 
as 


= ^g pv {n d v¥+ ¥Yvd„¥-d^¥Y v ¥-d v ¥Yn¥) 

- l - A g pv ¥(Yn r v + r v y M + 7 vF m +r M y v ) y/-S p (2 KF k -F{K)) . (2.13) 

Finally, exploiting the explicit form of spin connection (12.91) after some manipulations one finds 
the following non-trivial components of the energy-momentum tensor of the spinor field 


7^0 

7 0 

T\ 

t 3 

Ti 


m sp S + F(K), 

r 2 2 = r 3 3 = F(K) - 2 KFk, 


L o 1 Yd 1 _dA 

4 a i \a 3 a i / 
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A 2 , 
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(2.14a) 

(2.14b) 

(2.14c) 

(2.14d) 
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So the complete set of Einstein equations for a BI metric should be 


— + — + —— = k(F{K) - 2KF k ), 

a 3 a i <23 a\ 

2— + = k(F(K) — 2KF k ), 

a i a\ 

-j + 2—— = K(m sp S + F(K)), 
a\ a 3 ai 



(2.15a) 

(2.15b) 

(2.15c) 

(2.15d) 

(2.15e) 


Before solving the Einstein equations let us first write the equations for the bilinear spinor 
forms. Recalling that there are 16 bilinear spinor forms, namely, S = xfnj/, P = fy/y^y/. = 
y/yVy/, A“ = xfry 5 y u y/. and (2“ v = y/rA ,v ' y/ are the scalar, pseudoscalar, vector, pseudovector 
and antisymmetric tensor, respectively, for the LRS BI metric one finds the following system of 
equations: 


S 0 + &Aq = 0, (2.16a) 

P 0 - <Mq = 0, (2.16b) 

A® + d>Po — f^S 0 = 0, (2.16c) 

Aq = 0 , (2.16d) 

Vq = 0 , (2.16e) 

4 + ^ + ^Qo=0, (2.16f) 

q 3°-4*3 = 0, (2.16g) 

Ql l -&vl = 0, (2.16h) 


where we denote So = SV. Pq = PV, Aq = A^V. Vq = V, Qq V — Q^ v V and = m sp + We 
also denote = 2 SFkK] and (# = IPFkKj , with = dF/dK, Kj = dK/dF and Kj = dK/dJ. 

Here we also introduce the volume scale 


V = afa^. 


(2.17) 
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III. SOLUTION TO THE FIELD EQUATIONS 

From (12.16al) - (12. 16hl) one finds the following relations: 


(So) 2 + (F 0 ) 2 + (A[i) 2 = C 1 = Const, 

(3.1a) 

Aq = C 2 = Const, 

(3.1b) 

Vq — C 3 — Const, 

(3.1c) 

(^ 0 ) 2 + (e 5 1 ) 2 + (v^) 2 -C 4 = Const. 

(3.Id) 


Let us now go back to the Einstein equations. The off-diagonal components of Einstein equa¬ 
tions (12. 15dl) and (12.15el) impose the following restrictions either on the components of the spinor 
field or on the metric functions: 


A 2 = 0, A 1 = 0, 


03 

a 3 



(3.2a) 

(3.2b) 


The restriction (13.2bl) leads to a 3 — qocii with qo being some constant. In this case the system 
can be described by a Friedmann-Robertson-Walker (FRW) model from the very beginning. Here 
we do not consider this case, which we will address in a later work, within the scope of a FRW 
model. 

We consider the case when the restriction is imposed on the components of the spinor field in 
detail. Subtraction of (12.15bl) from (12.15al) gives 


a 3 a\ a\ \a 3 a\ 


(3.3) 


that leads to [5] 


ai =£>iV 1/3 exp (Xi f a 3 = (l/Dj)V 1/3 exp f-2Xi f ^ 


(3.4) 


V J 1 - v ' "V J V. 

with D\ and X\ being the integration constants. Thus we see that the metric functions can be 
expressed in terms of V. 

The solutions to spinor field equation (I2.10al) in this case can be presented as [5] 


= ^=exp^-j JQdt'J, ys 3 ,4(t) = ^=e x P (i j^dt^, 


C 3 ,4 


(3.5) 
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with Ci, C 2 , C 3 , and C 4 being the integration constants and related to Vo as 

j Cj + ^2^2 ^ 3^3 ^4^4 — vo* 

Here we assumed that K — I, i.e., F = F(7). The reason for this choice is discussed later. 

Thus we see that the metric functions, the components of spinor field as well as the invariants 
constructed from metric functions and spinor fields are some inverse functions of V of some de¬ 
gree. Hence any space-time point where V = 0 is a singular point. So it is important to study the 
behavior of V, which we do in the next section. 


IV. RESULTS AND DISCUSSION 


In this section we discuss the results obtained in the previous section. In doing so, we pay 
special attention to the volume scale, V. 

Let us first see whether the model becomes asymptotically isotropic. It can be shown that for 
an expanding Universe, when V —» °° as t —» °°, the isotropization process of the Universe takes 


place. To prove that we exploit the isotropization condition proposed in [26|] 


Oj. 

a 


—> const. 

t —>00 


(4.1) 


Then by rescaling some of the coordinates, we can make aija —> 1, and the metric will become 
manifestly isotropic at large t. 

Taking into account that a — V 1 / 3 from (13.41) we find 


~ = Di exp Jy)->D U y = (l/D 2 l )exp(-2X l j -A 1/D 2 , (4.2) 

as V —» 00 . Recall that the isotropic FRW model has the same scale factor in all three directions 
(i.e., a\{t) = < 22(0 = ^ 3 ( 1 ) = a(t)). So for the LRS BI universe to evolve into a FRW one we 
should have D\ — 1. Moreover, the isotropic nature of the present Universe leads to the fact that 
|Xi| << 1, so that —y 0 for t < °° (for V(t) — t n with n > 1 the integral tends to zero 


as t —> 00 ). 

Our next step will be to define V. Combining the diagonal Einstein equations (I2.15al) - (I2.15cl) 
in a certain way for V we find [5] 


V = — (m sv S + 2(F(K)-KF k )) V. 


(4.3) 
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Now in order to solve (14.31) we have to know the relation between the spinor and the gravita¬ 
tional fields. Using the equations (12. 1 6al) and (I2.16bl) it can be show that 


K 


V 2 

v 0 

V 2 ' 


(4.4) 


Relation (14.41) holds only for massless spinor field if K takes one of the expressions {7. I + 
7 ,1 — 7}, while for K = I it holds both for massless and massive spinor fields. In the case of 
K — I + J one can write S = sin (Vq/V) and P = cos (Vq/V), whereas for K = 7 — 7 one can write 
S — cosh (Vq/V) and P = sinh (Vq/V). In what follows, we will consider the case for K — I, as in 
this case further setting spinor mass m sp = 0 we can revive the results for other cases. Assuming 


on account of S = Vq/V we find 



(4.5) 


V 


3k 

T 


m sp Uo + 2 £4( 1 — n k )V* lk V'- 2rlk 
k 


(4.6) 


with the solution in quadrature 


dV 


3k \m sp V 0 V + ZkhVo lk V 2 ( 1 ~ n ^} +C 


— t + to, 


(4.7) 


with C and to being some arbitrary constants. 

Thus we see that the metric functions, the components of the spinor field, as well as the in¬ 
variants constructed from metric functions and spinor fields are some inverse functions of V of 
some degree. Hence any space-time point where V — 0 is a singular point. So we consider that 
the initial value of V (0) is small but non-zero. As a result for the nonlinear term to prevail in (14.61) 
we should have n^ — n i : 1 — 2n\ < 0 (i.e., n\ > 1/2) whereas for an expanding Universe when 
V — * °° as t —y oo one should have = nj : 1 — 2/72 > 0 (i.e., 772 < 1/2). As is seen from (14.61) . 
n k = no : no = 1/2 leads to a term that can be added to the mass term. 

In this case we obtain 


V = &i(V), 
^i(^) = y _ 


(m sp + Ao) V 0 + 2^ (1 - «!)V 0 2fil V 1 - 2 " 1 + 2A 2 (1 - nV)V^ 11 V 1 ~ 2 ” 2 


(4.8) 
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Equation (14.81) allows the first integral 
V = <S> 2 (V ), 


d> 2 (y) = . 3k [(/n sp + Ao) VoV + hV^'V 2 ^) +A 2 y 0 2n2 y 2 ( 1 -"2) +c 


(4.9) 

(4.10) 


The solution to (14.81) can be written in quadrature as follows: 

r dv 


® 2(V) 


— t+tQ. 


(4.11) 


To solve (14.81) we should choose the problem parameters Vo, m sp , k, C, A*-, as well as the initial 
value of y (0) in such a way that does not lead to 


(m sp +Ao) y 0 y+A 1 y 0 2 " 1 y 2(1 ~ ni) +A 2 y 0 2 " 2 y 2(1 '" 2) +c < o. 


For simplicity let us set Uo = 1, m sp = 1, Co = 10, and K = 1. In line with our discussion 
earlier we consider n 0 = 1/2, n\ — 2, and n 2 — 0. In our case we set V (0) =0.5. We set Ao = 1, 
whereas Ai = ±1 and k 2 = ±1 were taken in different combinations. It was found that depending 
on the sign of A 2 the model gives principally different types of solutions, namely, in the case of 
positive A 2 we have an accelerated mode of expansion of the Universe, while for negative A 2 we 
have oscillatory solution. 

Defining the deceleration parameter 


from (14.81) and (14.91) we have 


yy _ yd>i(y) 
<%(v) ’ 


(4.12) 


3 K 

2 

(m sp +Ao) V 0 V + 2Ai(l -n 1 )V 0 2 ' , 1 y 2(1 -' ,l) +2A 2 (1 -n^V^V 2 ^-'^ 

3k 

(m sp + Ao) y 0 y+AiV^y 2 * 1 -"!) + A 2 y 2 , 72 y 2 ( 1 - ,! 2) +c 



(4.13) 


Taking into account that for an expanding Universe at large t the term V 1 1,12 prevails, for 
deceleration parameter we find 


lim q —»— (1— n 2 )<0, since n 2 <l/2. (4.14) 

V—>°° 

Thus we see that spinor field nonlinearity generates late time acceleration of the Universe. 

In Fig. [H and Fig. [2] we plotted the evolution of volume scale V for a positive and negative 
self-coupling constant A 2 , respectively. As one sees from Fig. Q] a positive A 2 gives rise to an 
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FIG. 1. Evolution of the Universe for a positive Xi 
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FIG. 2. Evolution of the Universe for a negative A 2 


accelerated mode of expansion, whereas Fig. [2] with negative A 2 shows the oscillatory mode of 
expansion. In Fig. [3]the deceleration parameter q is illustrated for a positive A 2 . 

Finally we study what happens to shear and anisotropic parameters in this case. In doing so, let 
us first rewrite the corresponding quantities. The expansion tf for LRS BI metric reads 




a\ a 3 

2— H—- 
a 1 «3 


V 
V ’ 


(4.15) 







Spinor field with polynomial nonlinearity in LRS Bianchi type-I space-time 


11 



FIG. 3. Plot of deceleration parameter q for a positive Xi 


whereas, from 


i 9 a i 1 1 (ci\ cii 

a \ —°2— -~ ^ ~-P 

a 1 3 3 V a\ 113 J 


«3 


1 


2 (ti\ ci 2 , 


03 = ^ --- 

(23 3 3 \<3] <33 


one finds the expression for shear 


cr 2 = 


sGH* 2 


1 ( a \ (33 


3 \ <3i (23 


The anisotropic parameter in this case has the form 


1 


a, "- 3 £U 1 


Hi 


3 H 2 




at 


03 


1 , 


where H — ^ (2^j- + ^ ] = Further from (13.41) we find that 


a! 3 V + V ’ a 3 3 V V 
Now inserting (14.191) into (14.161) - (14.181) we finally find 


1 9 X \ 3 

>i = <yt = v , °3 


Xi 
V ’ 


y 2 

(j 2 = 3—i- 
y 2 ’ 


(4.16a) 

(4.16b) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 
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X 1 2 

A m = 18^. (4.22) 

As was shown in (14.81) and (14.91) in the case of a positive hi that generates late time acceleration, 
both V and V become large with the expansion of the Universe leading to <7/ —> 0, o 2 —> 0 and 
A m —> 0. This corresponds to our earlier conclusion regarding isotropization. As far as negative hi 
is concerned, in this case the model gives rise to an oscillatory mode of expansion. In this case we 
have both local minima and maxima. The maximum (minimum) value of volume scale V depends 
on the parameters and the initial condition and may be as large as possible. Hence in case of a 
negative A 2 though it is possible to attain a solution such that <7- |y=y max —> 0 and o 2 \v=v mm —> 0, 
but at the same time we have A m |y = v max —y °° because at any space-time point where V — V ma x(min) 
we have ^ |v=v maxfmmj = 0. This means that at any space-time point where evolution changes its 
direction (expansion to contraction and vice versa) the Universe becomes highly anisotropic. 


V. CONCLUSION 


Within the scope of the LRS BI cosmological model we studied the role of the spinor field in 
the evolution of the Universe. The reason for considering the LRS BI model lies in the fact that 
in the case of a full BI model the non-diagonal components of the energy-momentum tensor of 
the spinor fields imposes severe restrictions on the components of the spinor field, resulting in 


vanishing scalar S = y/y/ and pseudoscalar P — iy/ y 5 t/r [21 


22]. As a result both the mass term 


and nonlinear term in the Lagrangian disappear. But, as was shown here, in the case of an LRS BI 
cosmological model, neither the mass term nor the nonlinear term vanish. Moreover, unlike the 
Bianchi type-VI model the present model leads to asymptotic isotropization. It is was also found 
that, depending on the sign of the self-coupling constant, the model allows either the accelerated 
mode of expansion or the oscillatory mode of evolution. 
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